Compact quantum electrodynamics in 2 + 1 dimensions and spinon deconfinement: a 

renormalization group analysis 
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We discuss compact (2 + 1) -dimensional Maxwell electrodynamics coupled to fermionic matter 
with N replica. For large enough N, the latter corresponds to an effective theory for the nearest 
neighbor SU(N) Heisenberg antiferromagnet, in which the fermions represent solitonic excitations 
known as spinous. Here we show that the spinous are deconfined for N > Nc — 36, thus leading to 
an insulating state known as spin liquid. A previous analysis considerably underestimated the value 
of Nc. We show further that for 20 < A'' < 36 there can be either a confined or a deconfined phase, 
depending on the instanton density. For TV < 20 only the confined phase exist. For the physically 
relevant value N = 2 we argue that no paramagnetic phase can emerge, since chiral symmetry 
breaking would disrupt it. In such a case a spin liquid or any other nontrivial paramagnetic state 
(for instance, a valence-bond solid) is only possible if doping or frustrating interactions are included. 

PACS numbers: ll.lO.Kk, 71.10.Hf, 11.15.Ha 



I. INTRODUCTION 



Quantum electrodynamics in D-l-1 spacetime dimen- 
sions (QEDd i) with D = 1,2 are useful field-theoretic 
models in high-energy physics. Phenomena like chiral 
symmetry breaking and confinement are easier to under- 
stand in QEDi 1 and QED2_i than in QCD. The simplest, 
exactly solvable model of this type is spinor quantum 
electrodynamics in 1-1-1 spacetime dimensions (QEDi.i), 
the so-called Schwinger model, which exhibits both chi- 
ral symmetry breaking and confinement Another ex- 
ample, relevant to the present paper, is the (2 -I- 1)- 
dimensional spinor quantum electrodynamics, QED2,i, 
in the form introduced by Pisarski some time ago.^ This 
model is known to exhibit spontaneous chiral symmetry 
breakingi^ An interesting aspect of this model is its ap- 
plicability to condensed matter physics where it appears 
in different contexts, especially in the study of high- Tc su- 
perconductors and Mott insulators i^i^iiiSii^ In the study of 
Mott insulators, it arises as an effective theory of the so- 
called spin liquids, which are Mott insulators without any 
broken symmetry. In this case the Dirac fermions repre- 
sent the so-called spinous, soliton-like excitations carry- 
ing spin degrees of freedom but no charge. A good name 
for this type of QED is quantum spinodynamics (QSD), 
since it is actually a quantum field theory of spinous. 
The theory can be derived for Mott insulators, and it is 
found that the abelian gauge field coupling to the spinons 
is compact. This follows immediately by accounting for 
the fluctuations around mean-field theories of resonating 
valence bonds (RVB)iS states which have a local U(l) 
gauge freedom in which the phase angle is defined only 
modulo 2TT. These mean-field theories are derived from 
the strong-coupling limit of the Hubbard model, the so- 
called Heisenberg-Hubbard model: 



where Si are spin-1/2 operators formed from Fermi fields 
fil3 subjected to the local constraint fj^fia = 1 as 
Si — (1/2) f^^cTapfip, where cr are the Pauli matrices 
(T = {(Ti, (72, CT3) (throughout this paper summation over 
repeated greek indices is implied). The sum in ([T]) runs 
over nearest neighbor pairs of sites in a square lattice, 
and J is related to the original parameters of hopping 
and interaction energies t and U of the Hubbard model by 
J = At'^/U. Different composite-field theories are known 
to represent the same quantum states of the model. This 
follows from the fact that the local gauge symmetry of the 
Heisenberg-Hubbard model is actually SU(2)iii Thus, 
the composite link fields 



and 



Xij 



(2) 



(3) 



H 



(1) 



where are nearest neighbors, which are obtained 

from different Hubbard-Stratonovich decouplings of the 
Heisenberg-Hubbard model, describe the same physics 
when associated with the most stable ground state of 
the corresponding mean-field theory. They are connected 
through a SU(2) gauge transformation^ 

The phase fluctuations of either composite fleld fea- 
tures a lattice gauge fleld Aij. These are obviously 
of the compact U(l) type. Both fields together trans- 
form into each other by local SU(2) transformations, and 
theories utilizing this symmetry have been studied in 
the pastr^ii^ with some advantages over compact U(l) 
theories, especially if one is interested in studying the 
phase structure of high- Tc superconductors J^ii^ii^ A lo- 
cal SU(2) gauge theory also emerges in the study of frus- 
trated Heisenberg antiferromagnets i^^'^^ The existence of 
the larger SU(2) symmetry ensures the compactness of 
the Abelian theories even in the continuum limit, since 
the U(l) group is a subgroup of SU(2)ji^ 
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A controlled study of U(l) spin liquids was initiated 
some time ago by Affleck and Marston^ using the com- 
posite field ([Sj. To have a small expansion parameter, 
they generalized the (7/060? symmetry SU(2) to SU(N) and 
solved the model in the large N limit. This was done in 
the so-called self-conjugate representation of SU(N),— 12£ 
where the spin operators are given by Si^ap = f la fill ~ 
^ap/'i and fulfill the local constraint fj^fia = N/2. These 
operators have zero trace: Ti[Si^ap) = 0. If the partition 
function is represented as a functional integral over the 
action associated with the Hamiltonian ([1]) we may cal- 
culate a mean-field approximation from the saddle point 
approximation. The result which preserves all the lattice 
symmetries is the so-called 7r-flux phase, ^ whose spec- 
trum of elementary excitations is given by 



= 2|xo|yCOs2fc:r+COs2fcy, (4) 

where |xol is the mean- field amplitude of Xij- The excita- 
tions around the Fermi points ±(7r/2,7r/2) can be at low 
energies represented as four-component Dirac fermions. 
The phases of Xij fluctuate strongly and form a link gauge 
field Aij . It is not difficult to show that the effective low- 
energy Lagrangian in imaginary time has the formal 

1 ^ 

a=l 

where F^^, — d^A^j — d^A^ and we have used the usual 
Feynman slash notation = ^^a^j,, with 7^ being 4x47 
matrices. In the above continuum notation the compact- 
ness of the gauge field is not apparent. As it stands, the 
above Lagrangian is just the above-mentioned massless 
QED2,i, which is a well studied model. We shall reserve 
the abreviation QSD for the compact version of QED2,i, 
whose field theory will be discussed at lenght in Sect. III. 

In QSD the spinons play a role similar to quarks in 
QCD. Indeed, compactness of the gauge field leads to 
spinon confinement if N is not large enough. 
Note, however, an important difference. In QCD the 
gluon is introduced ad hoc to generate the coupling be- 
tween the quarks. The "gluon" in QSD, on the other 
hand, has a clear origin: it is spontaneously generated 
by the phase fluctuations of the link field ([21), which in 
turn is a composite field made of lattice fermions. This 
unique feature of QSD led Wen^^ to propose a similar 
mechanism in 3+1 dimensions to explain the origin of 
gauge bosons. 

Whether spinons in QSD deconfine or not was for 
some time matter of controversy<22i2SiiSii2^i22,!^ The con- 
troversy seems now to be solved at least at a qual- 
itative level. The result is that spinons deconfine for 
large enough A^,^^ and this guarantees the stability of the 
large- A'^ spin liquid. The aim of this paper is to further 
improve our understanding of spinon deconfinement, in 
particular by a more quantitative analysis based on the 
renormalization group (RG). We start in Sect. II with a 
discussion of the model in the absence of matter, which 



just corresponds to the compact version of Maxwell elec- 
trodynamics in 2-|- 1 dimensions (QED2,i).— The equiv- 
alence of this model with a three-dimensional Coulomb 
gas of instantons allows us to make a relatively simple 
field theoretical analysis, based on the equivalence of a 
a Coulomb gas with the sine-Gordon theory. In Sect. II 
we compute the one-loop effective potential of the model 
in terms of a scalar field 1^9(2;) whose correlation func- 
tion {ip{x)ip{x')) gives a direct measure of the interaction 
between instantons in three-dimensional spacetime. The 
calculation of the effective potential leads to an estimate 
of the instanton mass (which is different from the mass 
of the scalar field). We also derive the RG equations for 
the model. The reulting effective potential and the RG 
equations confirm in a physically appealing way the well- 
known result^^''^^ that in absence of matter the (2-1-1)- 
dimensional compact gauge theory permanently confine 
test charges. This result contrasts with the one in 3-1-1 
dimensions, where a deconfined phase has been proven 
to exist 

Section III is where we start with the analysis of com- 
pact QED2^i or QSD. Here is very difficult to write a 
field theory in terms of the original variables. It will 
be necessary, just like in Sect. II, to work with a field 
theory describing the dynamics of the instantons. We 
still have a gas of instantons, but it is no longer of the 
Coulomb type, since the interaction between the instan- 
tons is modified by the vacuum polarization. The RG 
equations are obtained in essentially the same way as in 
Sect. II, except that the effects of the vacuum polar- 
ization are taken into account. These change the phase 
structure of the theory in an essential way. By consider- 
ing a theory with A^ replica of Dirac fermions, we show 
that spinons deconfine for A^ > A'c = 36. This is a consid- 
erable numerical change with respect with our previous 
calculation of Ncr^ where a value smaller by a factor tt^ 
was found due to a wrong counting of tt factors in the 
RG calculation. Apart from this mistake, the analysis in 
Ref. [31I is correct, in the sense that a critical value of 
N is predicted above which deconfinement occurs. Here 
we elaborate further the nature of the confined phase. In 
particular we show that confinement is the only possibil- 
ity for A^ < 20 while deconfinement certainly occurs only 
for A^ > 36. In the region 20 < A" < 36 either phase is 
possible, depending on the instanton density. The uni- 
versal properties of QSD are encoded in the coefficient of 
the —1/R contribution to the interspinon potential. In 
the deconfined phase this coefficient is given by the fixed 
point in the renormalization flow derived in the 1 /N ex- 
pansion. However, in the confined phase the 1/A^ expan- 
sion does not hold and an expansion in the fugacity is 
made instead. Spinon deconfinement implies the stabil- 
ity of the spin liquid. The experimentally relevant case 
is the N — 2 one, which corresponds to the regime where 
the spin liquid is unstable, since the spinons are confined 
in this case. The nature of the confined phase for N — 2 
is discussed in Sect. IV, where the important role played 
by chiral symmetry breaking is emphasized. 



3 



II. COMPACT MAXWELL 
ELECTRODYNAMICS IN 2+1 DIMENSIONS 

The (2+l)-diinensional compact electrodynamics stud- 
ied by Polyakovfi^i abreviated here as QED2 i , is actually 
Maxwell theory in 2 + 1 dimensions in which the gauge 
group U(l) is made compact. The model was originally 
motivated by the spontaneous symmetry breaking pat- 
tern of the Georgi-Glashow model in 2 + 1 dimensions. 
After spontaneous symmetry breaking of the SU(2) group 
one is left with a residual U(l) group which is compact, 
since it is a subgroup of SU(2). Another place where this 
theory arises naturally is in lattice gauge theoryi^ 

It is well known that QED2 i is equivalent via dual- 
ity to a three-dimensional Coulomb gas of instantons,^- 
whose energy interaction is given by 



_ 2^^ 



(6) 



where cq is the bare gauge coupling and qi — ±q with 
q G N are the instanton charges. 

In two dimensions the Coulomb gas corresponds to 
a theory dual to the two-dimensional XY model^i^ 
or, equivalently, a two-dimensional classical superfluid.— 
The Coulomb interaction between the charges is in this 
case proportional to ln|xi — Xj\. In the context of two- 
dimensional superfluids, the charges in the Coulomb gas 
are interpreted as vortices in two dimensions. The two- 
dimensional Coulomb gas is known to undergo a vortex- 
antivortex pair unhiding phase transition, the cellebrated 
Kosterlitz-Thouless (KT) phase transitioni^ In the lan- 
guage of the Coulomb gas we have a low-temperature 
dielectric phase separated from the high-temperature 
"metallic" plasma phase by a phase transition without 
breaking the U(l) symmetry of the original superfluid 
or XY system, in agreement with the Mermin- Wagner 
theorem^S. In three dimensions, on the other hand, no 
phase transition occurs in the Coulomb gas and the sys- 
tem remains in the plasma phase. In this case the Debye- 
Hiickel (DH) mean-field theory is essentially correct and 
the interaction is screened in such a way that the ex- 
citations are always gapped. This can be conveniently 
expressed in field theory language by means of the sine- 
Gordon representation of the Coulomb gas. In the con- 
text of QED2 I the corresponding sine-Gordon theory 
reads 



1_ ,2 /^27r 
CsG = ■^(Ofj.y^) - 2zocos [—^ 



(7) 



The above Lagrangian corresponds to the field theory 
model dual to QEDj i- The parameter zq is the bare 
fugacity of the Coulomb gas. The DH theory amounts 
to a Gaussian approximation to the above Lagrangian. 
This leads to a mass Mq = ^-k ^f2z^ j et^^ for the scalar field 
This behavior of the dual model implies a correspond- 
ing mass gap in the magnetic field correlation functiouii^ 





27r2 


'^0 




3/2 




^ ~3~ 

















Fluctuation corrections to the DH approximation essen- 
tially do not change this result. To see this, let us com- 
pute the one-loop effective potential. This can be easily 
obtained by standard methods. At one-loop order it is 
more easily obtained by writing — (p 5lp, where (p is 
a constant background field while dp represents a small 
fiuctuation around it. By integrating out the Gaussian 
fiuctuations, the one-loop effective potential is obtained: 

Vcff(<^) -2zocos 

Veo / •i Lf^o V^o /J 

(8) 

where a counterterm proportional to cos(27r(/j/eo) 
was used to trivially subtract the contribution 
(2zoA/eo) cos(27r(/3/eo), with A being an ultraviolet 
cutoff. As with Eq. ([7]), the obtained effective potential 
implies a degenerate vacuum at p>n = neo, n e Z, whose 
energy density is given by Eq = — 2zo— (27r^/3)(zo/eQ)'^/^, 
which is always negative, just as the energy of the vac- 
uum without the quantum corrections. From this we 
conclude that the instantons of this theory are always 
massive. 

In order to better understand the meaning of this state- 
ment, let us compare the above effective potential with 
the one of a (l-l-l)-dimensional sine-Gordon theory. In 
this case the corresponding Euclidean theory is the dual 
field theory of a two-dimensional superfiuid,''^ which is 
known to undergo a KT phase transitioni^ The one-loop 
effective potential reads>^ 



-ttKz cos 



(2ttVKp ] In 



(2T:VK(pj 



2zcos I 2ttVKp 



(9) 



where K is the superfluid stiffness and A an ultraviolet 
cutoff. In order to remove the cutoff we add a coun- 
terterm proportional to cos {^.Ti^/Ki^ and impose the 

renormalization condition V^^{Q) = STr^ifz = Af^, which 
just fixes the renormalized mass of the scalar field in such 
a way as to have the same form as the one obtained from 
the DH theory. The result is 



Vcs{p) = -z{2 - ttK) cos (2TTVKipj 
-ttKz cos {2tt\/K(^ In cos {2tt\/Kp 



(10) 



The ground state energy density is now Eq — z{'kK — 2). 
This changes sign at Kc = 2/ti, precisely at the critical 
KT stiffness valuei^ In this case we can use the results 
of Ref. |43 to obtain the soliton mass as 



2 2z 
Msoi = -\ —{2-^K). 



(11) 



For K > Kc, which corresponds to low temperatures in 
the KT theory, the soliton mass vanishes while the energy 
becomes positive. For K < Kc, on the other hand, we 
can write 



sol 



8^0 



(12) 
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Note that while the sohton mass vanishes at Kc, the 
scalar field mass does not. This is very important, since 
in the (2+l)-dimensional case the mass of the scalar field 
also docs not vanish for any e^, but there the same is 
true of the mass of instanton excitations. By analogy to 
Eq. (fT2|) . we can infer that the instanton mass in CMT3 
should be given by 



Mnst OC - — 



eo 



IT- Zq 

^0 + — 1 — 

'0 



3/2' 



(13) 



which never vanishes in contrast to Eq. (I12p . 

The above results are made more transparent when we 
consider the RG equations of QED3, which are equivalent 
to the RG equations of the three-dimensional Coulomb 
gas. To lowest order this may be achieved by means 
of a scale-dependent DH approximation, in the same 
spirit of the analysis of the two-dimensional Coulomb 
gas made by Young4i The RG equations for the three- 
dimensional Coulomb gas were originally derived by 
Kosterlitz,"*^ using the so-called poor-man scaling ap- 
proach. As shall we see, for our purpose it is better to use 
Young's approach, "^^^ which we generalize to the higher- 
dimensional case. The derivation is given in Appendix A 
for the case of a d-dimensional Coulomb gas. Again, it 
will be useful to compare the RG equations for the KT 
case with those for QEDj i- Setting c? = 2 in Eqs. (|A17[) 
and (jAlSp of Appendix A, we obtain the well known RG 
equations for a two-dimensional superfluid where k = K: 



dK- 



dl 



dy 
dl 



TTK)y. 



(14) 



(15) 



From these equations we can see better why our approach 
of Kc from above drives the phase transition. In this case 
the sign on the right-hand-side of Eq. (fT5|l is negative, 
such that the system will flow to a regime of zero fugacity, 
leading to a line of fixed points. At this line the mass of 
the scalar field vanishes. If we define a dimensionless 
soliton mass through TOsoI = -^^soi/A, we obtain 



dm~ 



dl 



2(2/7r) 
K 



■[{2-TTK'f 



2y^K] 



(16) 



The fixed point of Eq. (I16p corresponds precisely to the 
KT critical point. In this equation it should be under- 
stood that for y — and K > Kc its right-hand side 
vanishes, since Mgoi is zero for all 2; > and K > Kc- 

Now let us see what happens in QED3 . In this case we 
set Kq = a = I/cq in the equations of Appendix A, where 
a is a short-distance cutoff, and define the dimensionless 
gauge coupling f = K^/ K = 1 / k = e'^ / Cq, such that 



/(/) = eQae'e(ae') = e'e(ae'). 



(17) 



with I = ln(r/a) being a logarithmic length scale. Thus, 
f{l)e~^ is just given by the screening constant ("dielec- 
tric" constant) of the Coulomb gas of instantons. In this 



way we obtain from Eqs. (lAlTp and (jAlSp the RG equa- 
tions for the QED3: 



dl 
dy 
dl 



vr 



y, 



(18) 
(19) 



Note that the coefficient of the term 1// in Eq. (|19p is 
different from the one obtained by us in Ref. i31i . where 
factors of tt were overcounted. It is easy to see that, in 
contrast with Eqs. ([Tl)) and (fT5)) . the RG flow pattern 
does not exhibit any fixed point, giving a further con- 
firmation that the three-dimensional sine-Gordon model 
does not undergo any phase transition. In other words, 
the excitations are always gapped. 



III. COMPACT QED IN 2+1 DIMENSIONS 
(QUANTUM SPINODYNAMICS) 

When fermionic matter is present, the energy of the 
instanton gas ^ is changed due to the vacuum polar- 
ization. The effective bare gauge coupling is modified to 
e^(p) = Z^(p)eg, where Za{p) is the wave function renor- 
malization of the gauge field in the noncompact theory. 
We have 



Za{p) = 



1 



i + n(p)' 



(20) 



where n(p) is the vacuum polarization. Thus, the energy 
of the instanton gas becomes 



where 



in' 



d^p 



(27r)3 p2 



-[i + n(p)], 



(21) 



(22) 



is the interaction between two instantons of opposite 
charge. 

Let us emphasize the dual nature of the above potential 
by writing the static semiclassical interspinon potential 
between two opposite test spinon charges: 



ViR) 



-el 



d^p 



(2^)2p2[i + n(p)]' 



(23) 



where R = |x|. We see that there is a manifest dual- 
ity between "electric" and instanton charges even after 
the vacuum polarization is included. Indeed, from Eq. 
((22)) the effective squared instanton charge is given by 
e^(p) = 47r2/e2(p), thus verifying the Dirac duality rela- 
tion e{p)e{p) = 2tt between the effective charges. 

At one- loop order we have Tl{p) = Ncg/Sp, implying 
that at large distances the potential becomes 



V{R) = - 



ttNR 



(24) 
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instead of having the logarithmic behavior expected clas- 
sically. Thus, quantum fluctuations lead in two spatial 
dimensions to a large distance behavior of the poten- 
tial similar to the one of electrodynamics in three spa- 
tial dimensions. Interestingly, the above power is \/R 
for all dimensions d E (2, 4) and the coefficient in front 
of the factor is universal, ■^^ a situation similar to 
the one obtained from string models for the interquark 
potential!^ 

Since the energy of the instanton gas is changed by the 
presence of matter, the sine-Gordon theory ([7]) must be 
changed correspondinglyi22i22iS^ It is easy to see that the 
field theory of the instanton gas is now given by 



C 



-^{d^,if )ZA{\/ ~d'^){d^ip) - 2zo cos(p, 
on 



(25) 



spinous, i.e., the spinous are deconfined and we have a 
spin liquid. On the other hand, for < 20 the instantons 
are in the "metallic" phase, which corresponds by duality 
to confined spinous. As we shall see, it is not quite accu- 
rate to say that spinous deconfine for iV > 20. Actually 
it will be shown later that in the interval 20 < < 36 
both confinement and deconfinement phases may occur. 
A stable deconfined phase occurs only for > 36. 

The arguments of the preceding paragraph are bet- 
ter understood by a RG analysis. We shall derive the 
RG equations using the method of Appendix A, with 
the Uo{r) given in Eq. (^5)1 replacing the one in Eq. 
(|Aip . The dielectric constant e(r) of the instanton gas 
as a function of the length scale r has the same form as 
in Appendix A, except that the renormalized interaction 
between the instantons is in this case given by 



where the symbol ZA{v~d^) has an operator mean- 
ing such that in momentum space the gradient energy 
reads simply {el/8TT^)ZA{p)p'^(p{p)(p{—p)- Thus, the La- 
grangian in Eq. (|25p corresponds to that in Eq. ^ 
with the bare squared charge Cg replaced by an effective 
one. Unlike in QED2 1 discussed in the previous Section, 
Eq. ([^5]) does not provide an exact dual representation of 
the theory, since its derivation is obtained in a harmonic 
approximation to the gauge field fluctuations. However, 
to the order to which our RG equations were calculated 
it provides an accurate dual field theory representation. 

After evaluating the momentum integral in (|22p using 
the short-distance cutoff a = 1/eg and an one- loop ap- 
proximation to n(p), we obtain the following interaction 
between two oppositely charged instantons: 



Uoir) 



TT 



An 



^ ^ W 2 N 



(26) 



where r = \x\. Let us compute the pair susceptibil- 
ity at large distances using the above potential. This 
is obtained from the r ^ 00 limit of Eq. (|A7p of Ap- 
pendix A for d — and inserting the Boltzmann factor 
n{r) « ZQe~'^"^'^\ The pair susceptibility obtained in this 
way gives a measure of (r^) and is given by 



Xo = lim Xoir) 



167r^z^ 
-34- 



dpp e 



(27) 



where the subscript in x indicates that the calculation 
is done with the potential UQ(r). By introducing the 
dimensionless variable u = Cqp/tt, we obtain 



16^8-^/4e 



Xo = 



H4r 



duu^-^/^e'/^. 



I/tt 



(28) 



The integral converges only for N > 20. This leads to 
£0 = 1 + 47rxQ > 1 signafing a "dielectric" phase for the 
instanton gas. Duality implies that the "dielectric" in- 
stanton phase corresponds to a "metallic" phase for the 



TT ds 



Uir) = Uia) + ^J^ + (29) 

which implies a renormalized dimensionless gauge cou- 
pling 



_ e'£(aeO 



(30) 



1 + —e- 

The first RG equation follows immediately from Eq. (I30p 



where we have defined 



2,,, _ 647r4 a6z2e6'-c^(«e') 



3 l + £e' 



(31) 



(32) 



The third term in Eq. ([31]) represents the correction to 
Eq. (|18p due to the presence of fermionic matter. For y ~ 
the instantons are suppressed and e = 1, such that Eq. 
((3T|) reduces to the /3 function of noncompact QED2.1. 
In such a situation we can rewrite the potential (|24p in 
terms of the fixed point = Att/N of the noncompact 
theory: 



ViR) = - 



(33) 



By differentiating Eq. (|32p with respect to Z, we obtain 



such that the correction to Eq. (|T9)) due to fermionic 
matter is also obtained. The flow of the screening con- 
stant of the instanton gas is given by 



de ey^ 
dl^~- 



This follows from 



de _ lOn^ 4a^ 4i_ui^^^i) 
dr ?, el 



(35) 



(36) 
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which is easily derived using Eqs. (jA6p and (|A7p of Ap- 
pendix A. 

The last term in Eq. (|34p is crucial for understanding 
QSD from a RG point of view. This term was overlooked 
bcforc^^'^'^-^'''-^^-^^ and led to an incomplete physical pic- 
ture of the phase structure of QSD. In order to emphasize 
the importance of this term, let us neglect it for a mo- 
ment. Considering /, = Att/N we find a line of fixed 
points for zero fugacity as N is variedj^i2^ From this 
it was conclude d^'^1^^ that QSD undergoes a KT-like de- 
confinement transition in three spacetime dimensions. In 
this deconfinement scenario the critical value of A'^ above 
which the spinous deconfine is Nc — 24. The reason why 
this term was overlooked before is that it was generally 
assumed that the logarithm term in Eq. (|26p dominates 
over the 3D Coulomb interaction at large distances. In 
other words, the large distance limit was being taken at a 
too early stage. We may argue that if one has a logarith- 
mic gas of instantons in 3D it should not be particularly 
surprising that a 3D KT-like transition emerges. How- 
ever, to really prove such a statement is far from being 
trivial and has led to recent controversies j^^J^i^ Monte 
Carlo simulations^ seem to support the scenario of a 
KT-like transition for the 3D logarithmic gas. However, 
it is now known that such a transition does not describe 
correctly the deconfinement transition in QSD. ^^'^^ Here 
we shall give a more complete explanation of the results 
obtained in Refs. [l^ and[3ll. 

In Ref. [l^ it was shown that the spinous deconfine for 
a large but finite N, but the critical value of N above 
which it happens was not calculated. This value was 
estimated in our paper'^^ to be Nc = From the 

more thorough analysis performed here we shall see that 
this estimate lies too low. One trivial reason for this is 
an overcounting of tt factors in the derivation of the RG 
equations, as was mentioned after Eqs. (US]) and 
But there are other aspects of the RG flow of QSD which 
were not considered in Refs. and and will now be 
included. 

From Eqs. ([?T|) . and ([55)) . we see that there is 

apparently for all A^ > 20 a fixed line at zero fugacity and 
/, = 47r£,/Af, 1 < e, < A^/20. Note that e > 1, since 
X > 0. However, it is clear that for vanishing fugacity 
we must have e, = 1, so that A = 20 and the fixed line 
actually collapses to a fixed point with ~ ir/b. It is 
reasonable to assume that the RG functions are given by 
expansions in the fugacity. For y — we have e{l) — 1 
for all I, such that Eq. (pij) reduces to the RG /3 function 
of noncompact QED3. Up to terms of order we can 
legitimately approximate /^/e ~ and yf/s ~ yf in 
Eqs. dniD and respectively. In this way Eqs. ^ 

and (IM)) decouple from Eq. (|35p and have the following 
fixed points at nonzero fugacity: 



/± = — (6±V36^ 
Af V 10 



(37) 
(38) 



We confirm once more the fixed point with zero fugacity 
at A^ = 20, y_ = and /_ — n/h. The above fixed points 
do not exist if A^ > 36, in which case only the fixed point 
of the noncompact theory is available, and the spinous 
are deconfined. This result confirms the analysis of Ref. 
[22l . where it was argued that spinous deconfine above a 
large but finite value of A^. Thus, for A^ > 36 the theory is 
well controlled by an expansion in 1/A^. For A^ e (20, 36], 
on the other hand, it seems to correspond to an expansion 
in the fugacity. In order to investigate this claim, let us 
consider Eq. (j30p in the limit Z — > cxd, in which case 
f ^ f± = 4:Tre±/N, where e± — limj^oo £(ae'). From 
Eqs. jST]) and dSSl), we obtain for 20 < A^ < 36 



and 



f± 



N 



1 



20^ + ;^^^ 



-yl 



(39) 



(40) 



The above equations clearly have the structure of an ex- 
pansion in y. This provides a further argument for setting 
e = 1 in Eqs. ([3T|) and (|34p. However, it should be noted 
that within the present approximation j/j. is not small in 
the entire interval 20 < A^ < 36, and perturbation the- 
ory may eventually breaks down for some values of A^. 
Most critical is the situation at the fixed point (/+, ?/+), 
since y+ gets small only for A^ close to 36. Note that also 
/+ can be large in this interval. Already at A^ = 20 we 
have /+ = TT, indicating that the behavior near the fixed 
point (/+, y+) should be considered with great care. The 
fixed point (/_,?/_), on the other hand, has much better 
asymptotics. Indeed, y_ is small for all values of A^ in 
the interval 20 < A^ < 36 up to A^ = 32, becoming larger 
than unity only above this value. 

For all 20 < A^ < 36 the fixed points {f±,y±) govern 
a confined regime such that the interspinon potential at 
these fixed points reads 



V±iR) = a±R~ 



f± 



-Vl{R) + 0{1/R^), (41) 



where a± is the string tension in the presence of fermionic 
matter. In the present approximation, the string tension 
has precisely the same form as in the QEDj 1,— except 
that the bare parameters are replaced by the renormal- 
ized ones. Thus, a = 2e^M/7r^ = 4e\/2z/7i', or in terms 
of dimensionless quantities, cr/eQ = A:\J2jy jix. The term 
Vh{R) is assumed to be given by a string model for the 
electric flux tube due to Liischer;^ which in d dimensions 
has the form 



y^R) 



(d- 2)7r 
24i? 



(42) 



Such a term has to be included in order to account for 
the fluctuations of the electric flux tube occuring even 
in absence of matter. The factor d — 2 represents the 
number of transverse modes of the string worldsheet. The 
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TABLE I: Summary of the phase structure of QSD. The third 
column gives the universal coefficient of 1/R. 



y 



N 


Phase 


— 1/R coefficient 


N <20 


confined 


7r/24 + /+/7r" 


20 < Af < 36 


confined/deconfined 


like iV < 20 or iV > 36 


iV > 36 


deconfined 


A/tvN 



coefficients of the l/R terms in Eq. (|41]l are universal 
and the string tensions a± reduce to the one obtained by 
Polyako\*i^ in the limit case where no matter fields are 
present. 

The schematic flow diagram for 20 < iV < 36 is shown 
in Fig. 1. In this case the system can be either in the con- 
finement or deconfinement phase, since the fixed points 
on each side of the dashed line in Fig.l are infrared stable. 
The fixed point (/_ ,y-) is infrared stable only along the 
dashed line in Fig. 1. Precisely for this reason it plays 
an important role in our discussion. The point is that 
vanishes for = 20, which is the value of A^ leading to 

= 1. This result gives a string tension that vanishes 
for all A^ < 20, thus leading to a collapse of the fixed point 
(/_,?/_) into the fixed point of the noncompact theory. 
Since the fixed point (/_ , y- ) is unstable in the direction 
leading to the fixed point and the one of the 

noncompact theory, the latter will no longer be stable 
and only the confinement phase governed by the infrared 
stable fixed point (/+, y+) will exist. Note that while the 
string tension a- vanishes continously as A^ approaches 
20 from above, this is not the case as A^ approaches 36 
from below. In fact, a±/eQ = '^\J1j±y±l-K has a finite 
value for A^ = 36. Since <j vanishes for A^ > 36, it follows 
that there is a universal jump in the string tension for 
AT = 36. 

It should be emphasized here that the physical inter- 
spinon potential corresponds to the one associated with 
the infrared stable fixed point, i.e., T^+(_R). Since the 
fixed point is unstable, the dashed line repre- 

sents a critical line separating confined and deconfined 
regimes. This role of a separatrix played by the dashed 
line of Fig. 1 is most clearly seen as N cross the value 
A^ = 20, in which case the line collapses along the axes / 
and y, with the fixed point becoming identical 

to the one of the noncompact theory. 



The phase of the system when 20 < A^ < 36 is de- 
termined by the size of the Debye-Hiickel parameter, 
Kb = nX^ = •\/2e^/(87r^Af ), where n is the instanton 
density and A_d the Debye length. Since for A^ > 36 
only the fixed point at zero fugacity exists, we have that 
cr+ = 0, i.e., deconfinement occurs for A^ > 36 and the 
large- A^ result is essentially correctj^ 

In Table I we summarize the phase structure of QSD. 



f 



FIG. 1: Schematic flow diagram for 20 < A'^ < 36 featuring 
the fixed points at nonzero fugacity given by Eqs. (I37|l and 



(EH). 



IV. DISCUSSION 



We have presented in this paper a description of the 
quantum spinon dynamics in a Heisenberg antiferromag- 
net using a field theory of A^ identical replica of fermions 
with a constraint. Only the limit N ^ oo can be treated 
exactly, the lower A^ results are obtained from different 
approximations. Since N = 2 corresponds to the phys- 
ically interesting case, one may wonder to which extent 
our results are physically relevant. This question is hard 
to answer. Originally, the study of spin liquids in 2-1-1 
dimensions was motivated by the properties of high- Tc 
cuprate superconductorsjiS, It was hoped that doping a 
Mott insulator provides a mechanism for superconduc- 
tivity in the cuprates. At zero doping the cuprates are 
antiferromagnetic insulators, i.e., the global SU(2) spin 
symmetry of the model is broken in the ground state. 
Spin liquids, on the other hand, have no broken symme- 
tries. Nevertheless, by doping a spin liquid we expect to 
obtain a superconducting stateiiii^ The idea was that 
doping may frustrate the magnetic state of the system 
and produce a spin liquid. With further doping the sys- 
tem would eventually become superconducting through 
a kind of hole condensation. Hence doping plays a role 
similar to frustration due to the triangular lattice in the 
organic compound k-(ET)2Cu2(CN)3, where a spin liq- 
uid phase was recently reported.>44 

Our analysis considers neither doping nor triangular 
lattices. Instead, the spin liquid emerges from the large- 
A^ limit, which provides the needed frustration to pro- 
duce a spin liquid. Our treatment at lower values of 
A^ was based on an expansion in the instanton fugac- 
ity. This treatment is essentially non-perturbative, and 
as such difficult to control. Nevertheless the stability of 
the spin liquid for large A^ seems to be established beyond 
doubt. 

Note, however, that the large-A^ analysis of the lat- 
tice model is made in a particular representation of the 
SU(N) group, by writing the spin operators in terms of 
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fermions.^'^^ If the spin operators are written in terms 
of bosons — which is an equally valid description of the 
problem — other representations are obtained. In such a 
case the Berry phase plays an essential role in determin- 
ing the phases of the system^^- There N is kept fixed and 
large, while a coupling constant g related to the exchange 
constant is varied. In the lower phase, i.e., for g < gc, 
where gc is the critical coupling, the SU(N) symmetry 
is broken and we have a Neel-like state. For g > gc, 
on the other hand, the paramagnetic phase is not a spin 
lic^uid. Instead, the so-called valence-bond solid (VBS) 
state emerges4s^ For N = 2 the VBS state can only oc- 
cur if frustrating interactions respecting the SU(2) sym- 
metry are included in the Hamiltonian. In such a sce- 
nario, where also an effective compact U(l) gauge the- 
ory is featured, the spinous are confined in both phases, 
but conjectured to be deconfined precisely at g — gc- 
This is the so called deconfined quantum criticality sce- 
nario and provides a new paradigm for quantum phase 
transitionsi^ One of the main characteristics of the this 
new type of quantum criticality is a large value for the rj 
exponent. In which systems this new paradigm actually 
holds is still under discussion. There is a recent numer- 
ical evidence^ supporting the scenario proposed in Ref. 
[50I . On the other hand, it was recently demonstrated us- 
ing large scale Monte Carlo simulations^^ that quantum 
antiferromagnets with easy-plane anisotropy do not ex- 
hibit deconfined quantum criticality, contrary claims in 
Ref. p. 

In QSD, our analysis certainly breaks down for N = 2. 
The description in terms of the generalized sine-Gordon 
theory ([25| no longer holds. The reason for this is that 
the derivation of the Lagrangian ([^5]) assumes that the 
fermions are masslesS j^^i^"^'^^ and it is known that in non- 
compact QED2.1 the fermion mass is dynamically gener- 
ated through chiral symmetry breaking (CSB) for small 
While the precise value of N below which CSB oc- 
curs is still a matter of debate, it seems to be more or 
less consensual that the chiral symmetry is broken for 
N = 2i^ Massive fermions modify the vacuum polar- 
ization in such a way that makes the derivation of an 
effective instanton Lagrangian difficult. Furthermore, in- 
dependent of the CSB occuring in noncompact QED2,i, 
confinement is likely to induce CSB in QSD already for 
N below and near 20, in a regime where the effective 
Lagrangian ([25|) can still be considered to be valid. In 
any case, CSB corresponds to spin density waves and 
this is precisely the state one would expect for undoped 
cuprates. The conclusion seems to be that no spin liq- 
uid state is possible for A'' = 2 in the undoped system, 
unless frustrating interactions are added. In this case it 
may well be possible that a VBS emerges instead a spin 
liquid. Moreover, the resulting paramagnetic state may 
depend on the nature of the frustration. Here it is worth 
to mention that a local SU(2) gauge theory of (fermionic) 
spinous for a frustrated Heisenberg antiferromagnet^^ ex- 
hibits a stable spin liquid phase. 

When doping is included the RG analysis becomes con- 



siderably more difficult due to the coupling of the gauge 
field with a non-relativistic complex scalar field.— In this 
case there are additional topological defects. These are 
vortex excitations coupled to the instantons. We are cur- 
rently investigating this situationi^ 



APPENDIX A: RENORMALIZATION GROUP 
EQUATIONS FOR THE d-DIMENSIONAL 
COULOMB GAS 

In order to make the paper self-contained we consider 
here the general d-dimensional Coulomb gas, whose RG 
equations were set up by Kosterlitz4^ They were origi- 
nally derived using the poor-man scaling approach. Here 
we employ the method due to Young,^ which is phys- 
ically appealing, since it amounts to applying a scale- 
dependent Debye-Hiickel argument which leads to the 
same results. Although Young applied the method to 
derive the RG equations associated to the Kosterlitz- 
Thouless (KT) phase transition, it can easily be gen- 
eralized to the d-dimensional case. We have done this 
previously^^ to derive the RG equations for anomalous 
Coulomb gases in d-dimensions. Here we concentrate on 
the ordinary d-dimensional Coulomb gas. 

The bare Coulomb interaction between two opposite 
charges of magnitude \/Kq is given by 



Uo{r) = -47T^KoV{r), 



where 



V{r) = 



, d 

r — 1 

4^-^/2 V 2 



2-d 



- 1 



(Al) 



(A2) 



In the above equation, a is a short-distance cutoff, which 
for d = 3 will be set to a — I/cq. From Eq. (jAip we 
obtain the bare electric field: 



Eo{r)^^AnMd)4^, 



where 



c(d) 



d-2 (d 



(A3) 



(A4) 



Next, in the spirit of the Debye-Hiickel theory, we in- 
troduce an effective medium via a scale-dependent dielec- 
tric constant e(r). This gives the renormalized electric 
field 



E{r) = -4TT^c{d) 



e(r)r 



(A5) 



The dielectric constant e{r) is expressed in terms of the 
the susceptibility x(r) as 



(A6) 
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where 



x{r) — Sd dss a{s)n{s) 



(A7) 



with Sd = 27r'^/^/r(d/2) being the surface of the unit 
sphere in d dimensions, and a(r) is the polarizability. 
For small separation of a dipole pair, it is given approxi- 
mately by 



a{r) 



(A8) 



The average number of dipole pairs is 

n(r) « z^e-^^'-', (A9) 

where zq is the bare fugacity and U{r) the rcnormalized 
potential obtained by integrating the rcnormalized elec- 
tric field (lA5l). 



£/(r) = U{a) + 'iTr^c{d)Ko 



ds 



la s'i-^eis)- 
The renormalized version of Kq is 



1 _ £(aeO 



Xd-2)l 



K{1) Ko 



(AlO) 



(All) 



where I — ln(r/a). Differentiating Eq. (jAlOp with re- 
spect to I, and using Eq. (|A11|) . yields 



dU ^fc{d) 



(A12) 



Next we differentiate Eq. (jAlip with respect to I to ob- 
tain 



dA_ ^ »7: b^z^^a ^2d/^c/(ae') ^ - 2)K-\ (A13) 



Here we define 



^2^^-j _ 87r~5'^Z0 ^2dl-U(ae') 



(A14) 



such that Eq. (jAlSp becomes 

dK-'^ 

~dr 



{d-2)K- 



(A15) 



From Eq. (|A14[) we derive the RG equation for the effec- 
tive fugacity: 



dz 



2TT^c(d)K' 



-,d-2 



(A16) 



ft is convenient to introduce the dimensionless quantities 
K EE a^^'^K and y ee a'^z to rewrite Eqs. (|A15|) and (|A16p 
as 



dn 



dl 



y^ + id-2)K-\ 



^ = [d-2nMd)^]y. 



(A17) 
(A18) 



For d = 2, the above RG equations govern the scaling 
behavior of the KT transition, while for d > 2 there is 
no fixed point, implying that the d-dimensional Coulomb 
gas is always in the metallic phase.**^ 
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